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Abstract 

A continuous theory which describes the process of doping in high- 
Tc superconductors of the YBCO-type is presented. It is shown that 
the introduction of dopants produces the creation of skyrmion topo- 
logical excitations on the antiferromagnetic background. The energy 
of these excitations is evaluated. Destruction of the Neel ground state 
occurs for a critical value 6c of the doping parameter for which the 
skyrmion energy vanishes. We evaluate this at zero temperature, ob- 
taining 5c = 0.39 lb 0.03 which agrees with the experimental value 
S^^P = 0.41 ± 0.02. 
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There is by now clear experimental evidence [jl], ^ Q that the pure high-Tc 
superconducting cuprates are well described by a quasi two-dimensional Heisenberg 
antiferromagnet on a quasi-square lattice, whose sites are occupied by Cu^^ magnetic 
ions. The continuous version of this |^ is the 0(3) nonlinear sigma field theory 
in three-dimensional spacetime and the ordered magnetization corresponds to the 
nonlinear sigma field ra", a = 1, 2, 3. This is subject to the constraint n°'n°' = p^, where 
^ = p'^—-^^=^ S being the spin, c the spin- wave velocity and a, the lattice spacing 
0. Using linear spin- wave theory results for the Heisenberg model, on the other 
hand, we find [Q] {^)lsw = 1.18{2S)^/2J where J is the Heisenberg antiferromagnetic 
coupling. These two relations are valid for large S. We assume the ratio of both holds 
for any S and, in particular, for S* = 1/2 we establish 

pI = 1-18^J (1) 

In this work we start from this continuous description of the pure system and propose 
a model for the doping process in YBCO-type materials, already in the continuous 
version. A picture emerges, then, that the doping process produces the formation of 
skyrmion topological defects on the antiferromagnetic background, thereby reducing 
the ordered magnetization. The skyrmion energy for an arbitrary doping is obtained 
through the evaluation of the quantum skyrmion operator correlation function and 
subsequent analysis of its large distance decay. The skyrmion energy works as an order 
parameter for antiferromagnetic order and this is destroyed when that vanishes. This 
allows us to obtain the critical value of the doping parameter at zero temperature. 

The nonlinear sigma field is conveniently expressed in the so called CP^ language 
in terms of a doublet of complex scalar fields Zi,i = 1,2 subject to the constraint 
zjzi = pg, namely 



= -4^tj^j (2) 
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where a" are Pauli matrices. In this language, the continuous field theory correspond- 
ing to the Heisenberg antiferromagnet is described by the lagrangian density in 2-1-1 
dimensions ^ 

Ccp^ = {D^Zi)\D'^Zi) (3) 



where = 9^ + iA^ and = ji{iz\d^Zi). 

The CP^-nonhnear sigma model possesses classical topologically nontrivial solu- 
tions called skyrmions. These bear one unit of the topological charge Q = J d?xJ^, 
where is the identically conserved topological current. In CP^ language = 
Q^j^p and we see that Q is nothing but the "magnetic" flux of the field A^. 
This means that, in this language, the skyrmions are vortices. In the nonlinear sigma 
language, they are defects in the ordered Neel state. A general method of vortex 
quantization ^ which allows the explicit construction of a vortex/skyrmion cre- 
ation operator therefore, can be applied to the skyrmions appearing in this theory. 
Before doing it (Eq. ([TB|)), let us consider the process of doping. In order to describe 
it we already start with the continuum language and associate to the electron/hole 
dopants a two-component Dirac field with two internal indices corresponding to the 
spin orientations: 'ipa{x),a =t, [. For Y Ba2Cu-iOQ^s the Fermi surface with dopants 
has an almost circular shape [0] and expanding around it we can use a relativistic 
approximation for the dispersion relation. This is clearly not valid for La2-sSrsCu04 
because of the shape of the Fermi surface 0. We propose the following lagrangian 
density for describing the dopants and their interaction with the backgroud lattice in 
YBa2Cu30(i+s 

where m* and vp are respectively the effective mass and Fermi velocity of dopants. 
This must be supplemented by the constraint zjzi = p^. The dopant dispersion 
relation corresponding to (H) is e{p) = \Jp'^Vp + {m*Vpy which agrees with the fact 
that the system is still an insulator even after the inclusion of some doping. The spin 
of the dopants is 5 = ipiBijipj. Notice that the ordered spin of the Cu^^ electrons 
is given by @ and their coupling to the ^^-field is described by iz\ 0^ ZiA^. This 
suggests the minimal coupling of the dopants to A^, introduced in @j. 

The doping parameter 5 is introduced by means of a constraint in the fermion 
integration measure, namely, 

(5) 
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where A'^ = 6 Jxl d^^5^{z — ^) for a dopant introduced at the point X = {X, T) and 
moving along the hne L. The 0*'*-component, for instance, is A^{z, t) = 6 6'^{z—X{t)). 
The above constraint can be implemented by integration over a vector Lagrange mul- 
tiplier field coupled as 



C\ — X, 



(6) 



In order to complete our description of doping, we must introduce their electro- 
magnetic interaction. The real electromagnetic interaction of particles of charge e and 
associated to a current = ipal'^'^Pa in 2+1 dimensions is governed by the lagrangian 
density [|ll[ 



EM 



-If. 



(7) 



which, therefore, must be included. 

Before collecting the pieces Cz,i),£x and Cek-u let us remember that the space of 
configurations of the CP^-nonlinear sigma model has inequivalent topological classes 
characterized by the Hopf invariant, Sh-, and corresponding to the nontrivial mapping 

— > S"^. When performing a functional integration over the Zj-fields, therefore, 
we must include in principle a weight exp[26'S'j^] where 5 is a real parameter to be 
determined. This implies the addition of a Hopf lagrangian density Ch, which in 
CP^-language is given by 

Ch = ^-^^'^''A.d^Ap (8) 

It has been shown that ^ = for the pure system however, in the presence of 
doping one has, in principle to add Ch well. The complete lagrangian density 
proposed for the description of YBa2Cu30Q+s is therefore 



^ — ^z,ip + -^^A + Cem + C 



-H 



(9) 



This must be supplemented by the CP^-constraint on the Zj-fields. 

We remark at this point that the theory described by (H) is only valid up to the 
threshold where the system is in the antiferromagnetic phase corresponding to a Neel 
ground state. Above this point there is no longer an ordered magnetic moment and 
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therefore the CP^-nonhnear sigma constraint is no longer justified. We do not ex- 
pect, therefore, that the superconducting transition could be explained by (|) without 
modifications. Observe also that in the region of validity of the presence of an 
insulating gap is clearly justified. 

We now functionally integrate over the fields Zi, ipa, and A^, in order to obtain 
an effective lagrangian for Afj,. The Zj-integral is easily done in the approximation 
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where = = constant {pl^ + pl^ = p1) yielding a piece —\^^lu[■zy\^^^ for 
the effective lagrangian density. The integration over the fermionic fields ipa is done 
by evaluating the fermionic determinant for an external field A = A^ + + eA^. In 
the small effective mass (m*) approximation this yields [|T^] 

]^e^'^PA,[T]d^Ap-\F,^mF^^ 

where F and S are kernels whose Fourier transform is F = + " ""Zl^ and S = 
sTP ~'~ "l.irik'^) • '^^^ remaining integrations over and A^ are quadratic and can be 
performed exactly thereby leading to the effective lagrangian density for the ^^-field: 

^^s[A^] = y-'^A^d^Ap - \:f^A^]^'''' - A/.^'^ + ^A/. [^P^" + AC^I K (10) 



where P^'^ = —OS^^'^ + d'^d'^, C^" = e^°"^da and 11 and A are kernels whose Fourier 

a+e'^/hc _ jni'vp/h i a _ _^ _ am'vp/ 
(A:2)3/2 fc4 aiiU iV — ^2 (A;2)3/2 



transform is 11 = ",'1.2 J Jl^ ~ ^"^ ^^^^ and A = — ^tt^w^- previous expressions. 



a,'y, P,a are numerical coefficients given by a = ^a^^^g , 7 = ^^S^iep^' ~ '"^'^ 



_ 327r^(24-7r^) 
^ - (7r2+16)2 • 



The field equation corresponding to (|T0]) is 



ee'^'^^d^Afs = A'^ + p2 



A^" + (—) d'^id^A") 



'111 



-□/ 

Taking the 0*^-component of it and applying to the skyrmion we see that the last 
term vanishes and 

ee'^diAj = A° or 68 = 6 6\z- X{t)) (12) 

where X{t) is the dopant position at a time t and B is the "magnetic" fiux or vorticity 
of A^. For the skyrmion B = 6'^{z — Xs(t)), where Xs{t) is the skyrmion position. 



We therefore conclude that the skyrmion position coincides with the dopant position 
at any time and 6 = 6. We conclude that upon doping, skyrmions are created on the 
ordered Neel background. This process reduces the ordered magnetic moment and 
eventually leads to the complete destruction of the antiferromagnetic order. A similar 
situation occurs in polyacetilene where the doping process produces the formation of 



solitons 0]. Observe that 6 vanishes for zero doping, in agreement with ||10| . 

Let us evaluate now the quantum skyrmion correlation function. Applying the 
formulation developped in p|, ^ to the system described by (0), we get 



< Kx)iJ,\y) >s= Z 



-1 



exp — Ss[J'fj,i> + B^jy] + Sqp 



(13) 



where Ss is the action corresponding to ([T^), Sqf is a. gauge fixing term and B 
2^ Ix L d'C.x^''^^'^ 6^ — ^) is an external field corresponding to a unit "magnetic" flux 
along L 1^, 0]. In order to put ( ]TD| ) in the form Cs[J-'^u] we wrote, up to gauge terms 
~^ Also in ([131) , we used 6 = 6 and identified the skyrmion and dopant 

positions: x = X, y = Y, according to (pT]), (|l^) and the discussion thereafter. 
Performing the quadratic integration over in (p!3D and neglecting the unphysical 
self-interaction LL-terms we get 



< lJ'{.x)ij){y) >-- 



exp 



1 

he 



6'^'ym*vpe 
8^r 



\x - y\ 



6\a + eyhe) 



In \x — y\ 



-i—[arg{x - y) + tt] 
/vr 



(14) 



From the large distance behavior of this we can extract the skyrmion energy for an 
arbitrary doping: 



r2 7^C 



(15) 
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In the last expression, we have written the Fermi momentum of the dopants as 
m*vp = where ap) is the lattice spacing for the dopants. The skyrmion en- 
ergy is an order parameter for the antiferromagnetic order and we conclude that this 
is destroyed for a critical doping parameter 



5c 



\ 'yhe 



(16) 
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The experimental values of the Heisenberg coupling and spin-wave velocity for 
YBa2Cu30Q+5 are, respectively, J = 100 ± 20 meV |, ||, and c = 1.0 ± 0.05 ^ 
0, 0. From J, we get pg through (|l]) . The only remaining input is the dopant lattice 
spacing. We take this to be the spacing between the oxygen ions in the CM02-planes, 
which for this compound, is ao = 2,68 A. Entering these data and the numerical 
value of 7 (after Eq. (|l^)) in (|^) we get 

5c = 0.39 ± 0.03 (17) 

This agrees with the experimental value of the critical doping at zero temperature 
namely 

= 0.41 ± 0.02 (18) 

We conclude that the continuous theory given by (0) provides a good descrip- 
tion of YBa2Cu30Q^s for ^ < ^c^ i-e., within the antiferromagnetic ordered region. 
The results also point towards a scenario in which holes are introduced in the p- 
orbitals of the oxygen ions in the Cu02 planes through doping. For a compound like 
La2-.sSrsCu04^, for instance, the relativistic approximation certainly is not valid and 
we do not expect a quantitative agreement with the experiment. Nevertheless, it is 
plausible that the mechanism of skyrmion formation by doping holds in general and 
is responsible for the destruction of the antiferromagnetic order in superconducting 
cuprates. The theory is presumably useful for 6 > 6c a.s well, provided some suitable 
modifications are introduced. We are presently investigating this possibility. 
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